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ABSTRACT 

O , We present predictions for the cosmic metal budget in various phases of baryons 

from redshift z — 6 — > 0, taken from a cosmological hydrodynamic simulation that 
includes a well-constrained model for enriched galactic outflows. We find that substan- 
tial amounts of metals are found in every baryonic phase at all epochs, with diffuse 
intcrgalactic gas dominating the metal budget at early epochs and stars and halo gas 
dominating at recent epochs. We provide a full accounting of metals in the context 
of the missing metals problem at z w 2.5, showing that ~ 40% of the metals are in 
OO ' galaxies, and the remainder is divided between diffuse IGM gas and shocked gas in 

\Q \ halos and filamentary structures. Comparisons with available observations of mctal- 

£N| . licity and metal mass fraction evolution show broad agreement. We predict stars have 

00 ' a mean metallicity of one-tenth solar already at z = 6, which increases slowly to one- 

, half solar today, while stars just forming today have typically solar metallicity. Our 

Hi column density- weighted mean metallicity (comparable to Damped Lya system 
metallicities) slowly increases from one-tenth to one-third solar from z = 6 — > 1, then 
C"| . falls to one-quarter solar at z — 0. The global mean metallicity of the universe tracks 

~ 50% higher than that of the diffuse phase down to z ~ 1, and by z — it has 
' \ a value around one-tenth solar. Metals move towards higher densities and tempera- 

tures with time, peaking around the mean cosmic density at z = 2 and an overdensity 
of 100 at z = 0. We study how carbon and oxygen ions trace the path of metals 
in phase space, and show that O III-O VII lines provide the most practical option for 
constraining intergalactic medium metals at z < 2. 
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1 INTRODUCTION 

Heavy elements are produced exclusively in stars, and 
then distributed throughout the Universe through su- 
pernovae and stellar winds. As such, they provide a 
unique tracer of the star formation and feedback pro- 
cesses that are central to understanding how galaxies 
form and evolve. Metals are present at the earliest ob- 
servable epochs: The farthest known quasar at z = 6.4 
shows substantial enrich ment despite being less than 1 Gyr 
old jBarth et al. 1120031) . and metals are seen in the in- 
tergalactic medium (IGM) at z > 
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lower redshifts, the gal axy mass-metallicity relation has been 
quantified out to ~ 2 iTremonti et al. II2004L ISavaglio et al. I 



1200.1 lErb et, al~ll200( ), intrac luster gas metallici ties have 
been constrained out to z ~ 1 iRosati et al. 1 l2004l) and the 
IGM is observed in metal lines nearly continuously from 
z ~ 6 -> 0. 

IPettini I 1^)99) first noted the interesting fact that the 
metals seen in Lyman break galaxies at z ~ 2 — 3 fall well 



short of the amount expected to be produced by the stars in 
those galaxies. This is known as the missing metals prob- 
lem. It implies that most metals must either reside out- 
side of galaxies, or be in some non-observable form in or 
around galaxies. Attempts to quantify m etals in the IGM 
iFerrara. Scannapieco fc Bergeron 1120051 . hereafter FSB05) 
or other galaxy populations such as s ub-millimeter galax- 
ies l|Bouche. Lehnert fc Pero ux ll20 06aT) or damp ed Lyman 
alpha (DLA) systems iProchaska et al. I I2006T) have not 
clearly been able to account f or the shortfall, though un- 
certain ties remain substantial dBouche. Lehnert fc Peroux I 
2006b, hereafter BLP06). Virtually every galaxy class and 
environment examined seems to contain a small but non- 
trivial amount of metals, so it appears that metal pollution 
has been curiously democratic. As the distribution of met- 
als in various environments and at various epochs becomes 
better measured, this provides important constraints into 
galactic feedback processes. 

A theoretical framework for understanding these wide 
ranging observations is provided by simulations of galaxy 
formation that include metal production and distribu- 
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tion m echanisms. The pioneering work of ICen fc Ostriker"! 

using cosmological hydrodynamic simulations high- 
lighted some basic trends arising in such models, such as 
the early enrichment of galaxies and the natural establish- 
ment of a metallicity gradient with cosmic overdensity. Using 
more advanced, higher-resolution cosm ological simulations, 
iThacker. Scannapieco. fc Davis I (I2002T) studied global metal 
enrichment from dwarf galaxies down to z — 4, and showed 
that it is possible to enrich the IGM at early times to a value 
broadly consistent with data while producing dwarf galaxies 
with enrichment levels similar to local dwarfs. 

An alternativ e ap proach was taken by 
ICalura fc Matteucci I l|2p04), who used detailed chemi- 
cal evolution models for various galaxy types together with 
observations of luminosity functions in those types to con- 
struct the evolution of metals in galaxies, intracluster gas, 
and IGM gas. While their models make detailed predictions 
for individual elements in various galaxy classes, they do 
not embed the predictions within an ab initio hierarchical 
structure formation framework, relying instead on passively 
evolving observed luminosity functions back in time. Still, 
this provides some interesting pred ictions that are not 
dissim ilar to those presented here. In ICalura fc Matteucci I 
( 2006) they extend their approach to study IGM enrich- 
ment, and suggest that winds with velocities exceeding 
1000 km/s are required to pollute the IGM to observed 
levels. 

In IQppenheimer fc Pavel (120061 hereafter OD06) we 
took such studies to the next level of sophistication, using 
cosmological simulations that incorporate enriched kinetic 
feedback in a full hierarchical structure formation setting. 
Our models, employing outflow scaling relations taken from 
local starburst data, showed remarkable success at repro- 
ducing observations of C IV enrichment in the IGM from 
z ~ 6 — + 2. Matching these data generally required highly 
mass loaded but low velocity winds from early galaxies, so as 
to eject a substantial metal mass into the IGM without over- 
heating it. Such scalings naturally arise in momentum-d riven 
wind scenarios iMurrav. Quatert. fc Thompson 112005)) that 
are also favored by local outflow obs ervations iMartin I 
120051: iRupke. Veilleux fc Sanders 1 120051) . providing an in- 
triguing connection between o u tflows at all epochs. In 
iDave. Finlator fc Oppenheimer I £2006) we show that sig- 
nificant early feedback is also necessary in order to sup- 
press reionization-epoch star formation. In Finlator, Dave 
& Oppenheimer (in preparation) we show that the result- 
ing galaxy mass-metallicity relation from our favored out- 
flow models agrees well with observations at z ~ 2. These 
successes suggest that we now have a plausible structure for- 
mation framework with which to study the global evolution 
of metals in the Universe. 

In this paper we use cosmological hydrodynamic simula- 
tions incorporating our most successful feedback model from 
OD06 (described in H2.H to study the evolution of metals 
in various baryonic phases in the Universe from z — 6 — » 
( H2.2I and provide quantitative answers for the missing met- 
als problem f H2.3l . We go on to compare the metallicities in 
those phases with observations f il2.4t . showing broad agree- 
ment. We make predictions for how metals evolve in cosmic 
phase space ( tj3.ll . and study how different oxygen and car- 
bon ions may be used to probe various regimes of phase 
space f H3.2I . We present our conclusions in jgl 



2 METAL EVOLUTION 
2.1 Simulations 

We run a cosmological hydrodynamic simula tion using 
Gadget-2 dSpringel fc Hernauist I l200l l2003ah with im- 
provements as described in OD06. Our run has 2 x 256 3 
particles in a randomly-generated volume of (32/i _1 Mpc) 3 , 
with an equivalent Plummer softening length of 2.5h~ kpc 
(comoving). Initial con ditions were generated using an 
lEisenstein fc Hul (Il999|) power spectrum , assuming a con- 
cordance cosmology Jspergel et al. 1 12003ft having S7 = 0.3, 
A = 0.7, n = 1, Ho = 70 kms" 1 Mpc" 1 , and a 8 = 0.9. 
The simulation was started at z — 159, in the linear regime, 
and evolved to z = 0. We include metal-line cooling by aug- 
m enting the pr i mordi al cooling rate calculated as described 
in iKatz et al. I <ll996f) with the difference in cooling rates 
between the primordial and metal-enriched collisional ion - 
ization equilibrium models of ISutherland fc Dopita I il993l) . 
Star forma tion is modeled in a multi-pha se manner as de- 
scribed in ISpringel fc Hernauist I (12003311 . occuring m re- 
gions sufficiently dense to undergo multi-phase collapse, 
with the threshold density depending on the metallicity 
through the cooling rate. 

We employ a momentum-driven wind model as de- 
scribed in OD06, in particular the vzw model, which pro- 
vides the overall best match to the C IV observations con- 
sidered. In our heuristic outflow model, a star forming 
particle has a probability to enter a wind that is pro- 
portional to the assumed mass loading factor, and if se- 
lected it is given a velocity kick. In the vzw prescrip- 
tion, the wind velocity is proportional to the galaxy ve- 
locity dispersion, which we estimate from the local value 
of the gravitational potential. The value of the velocity 
kick has some spread associat ed with an observed spread 
in th e the luminosity factor |R upke 1 Veilleux fc Sanders I 
2005); see IQppenheimer fc Dave I J2006h for full formula. 
As predicted in a momentum-drive n wind model (e.g. 
iMurrav. Quatert. fc Thompson 1120051) . we assume that the 
mass loading factor scales inversely with the velocity disper- 
sion. 

We identify galaxies using Spline Kernel Interpolative 
Denmax (SKID), and hal os using a spherical o verdensity 
algorithm, as described in lFinlator et al~l J2OO6) . Our sim- 
ulation galaxy mass resolution limit is M* > 9.91 x 1O 8 M0 
(i.e. 64 star particles); all predictions regarding galaxies are 
implicitly for galaxies above this stellar mass. 

Our simulations assume a constant metal yield of 0.02, 
hence metal production directly tracks star formation. Oth- 
ers prefer a value of 1 /42 to relate star fo rmation to metal 
mass production fe.g. iMadau et al~lll996t) . Most of the re- 
sults presented here will not depend sensitively on this 
choice, and in fact as long as this value is constant with 
time, the fractions of metals in various phases will be essen- 
tially unaffected (modulo small differences due to metal-line 
cooling). We assume a solar metallicity value of 2% metals 
by mass. 

While we focus on the vzw model here, the mzw model 
(the other momentum-driven wind model favored by OD06) 
produces similar results. We choose the 32/i _1 Mpc vzw run 
for this work because it has sufficient resolution to capture 
early star formation, while still having enough volume to be 
meaningfully representative at z = 0. On a practical note, it 
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Figure 1. Top: Metal mass fraction in various baryonic phases as 
a function of redshift: Stars (black solid), star-forming gas (blue 
dotted), halo gas (green dashed), shocked IGM (T > 10 4 - 5 A"; 
red long-dashed), and diffuse IGM (cyan dot-dashed). Middle: 
As above, for total mass fraction of baryons. Bottom: As above, 
for metallicity. The thick cyan line represents the mean mass- 
weighted metallicity of the universe. 



is the only model we have evolved to z — 0, owing to the sig- 
nificant computational expense involved. We note that cos- 
mic variance could still be significant, and the wind model 
is as yet not tightly constrained, so the exact values quoted 
here should be taken as illustrative rather than precise pre- 
dictions. Using a (16/i~ 1 Mpc) 3 volume vzw run, we have 
checked that resolution effects introduce only minor devia- 
tions from the numbers quoted herein from z = 6 — > 1.5. 

Although we are less confident of our predictions at 
2 < 2 because this particular feedback model has not been 
critically tested there, we include predictions anyway to il- 
lustrate some basic physical trends that are likely to be valid. 
As we will show later, the comparisons to data at z < 2 seem 
to indicate broad agreement. 



2.2 Metal Mass Fractions 

Figure (top panel) shows the evolution of the metal mass 
fraction in various baryonic phases: Stars, star-forming gas, 
halo gas, shocked IGM, and diffuse IGM. Middle panel shows 
the mass fraction in baryons for comparison. Star-forming 
gas corresponds to gas that is abov e the density threshold 
for multi-phase fragmentation in the lSpringel fc Hernauistl 



J2003al) recipe; this may be thought of as insterstellar 
medium (ISM) gas. Halo gas is defined as gas that is within 
identified halos but is not star-forming; this includes e.g. hot 
intracluster gas. Shocked and diffuse IGM gas encompass all 
gas outside of halos, a nd are separated b y a temperature 
threshold of 30,000 K JPave et al. Ill999l) . We choose this 
threshold as oppos ed to the canonical warm-hot interga lactic 
medium (WHIM: ICen fc Ostriker1ll999l:ll5ave et al. lEoOlh 
threshold of T > 10 5 K because in our models, substantial 
amounts of C IV absorption comes from gas right around 
10 5 K (OD06), whose temperature is set not by photoion- 
ization but by shock heating from outflows and/or gravi- 
tational collapse. Hence this gas rightfully belongs in the 
shocked IGM component. 

From z ~ 6 — * 1.5, diffuse gas is the largest reser- 
voir of metals; below this, stars take over with halo gas 
not far behind. It may be surprising that at early times 
the diffuse phase harbors most metals, given that metals 
are produced in much denser regions. This arises because 
early blowout is quite effective in our models, which is re- 
quired to both enrich the IGM early on and suppress early 
star forma tion to match observed high-redshift luminosity 
functions iDave. Finlator fc Oppenheimerl l2006l. The dif- 
fuse IGM contains virtually all baryons at high redshift 
as is generically seen in hydrodynamic simulations (e.g. 
iDave et al. Il200ll) . then drops steadily to around 30% mass 
fraction by today. The drop in metal fraction tracks the drop 
in total baryon fraction in this phase. 

As galaxies grow and outflows become less effective at 
polluting the diffuse IGM in our momentum-driven wind sce- 
nario, the metal mass fraction that remains locked in stars 
increases. By z ~ 2, stars contain 20% of all metals, and by 
today they hold over 40%. This value is consistent within 
uncertainties to the value o f 52% predicted in the model of 
ICalura fc M attcuccil i2004l) . Note that the stellar baryon 
fraction at z — is 6.8 %, which is in agr eement with obser- 
vational estimates (e.g. iBell et al. II2003T) . 

The amount of metals in shocked IGM gas remains 
fairly stable at 15 — 20% at all epochs. The mass fraction of 
baryons in this phase grows almost linearly with redshift up 
to 45% today. Hence while the shocked IGM contains a sig- 
nificant fraction of all present-day baryons, it does not hold 
a commensurately large reservoir of metals. Meanwhile, the 
growth of larger potential wells and reaccretion of IGM met- 
als into halos results in a gradual increase of metal mass in 
halo gas, up to around one-quarter of all metals today. Halo 
gas is generally over-enriched relative to its mass fraction. 

Star-forming gas shows a mild increase in metal content 
at early times, but after z = 3 it drops steadily to a very 
small number at the present epoch. This traces its baryon 
fraction evolution, which is difficult to see on this plot be- 
cause it is quite small at all epochs. This shows that galaxies 
are in general becoming more gas-poor with time, as inflow 
shuts off because of cosmological acceleration combined with 
the growth of large potential wells. 

In summary, a plausible model of cosmic enrichment 
that broadly matches key IGM and galaxy metallicity mea- 
sures predicts that metals are present in a range of baryonic 
phases from the earliest epochs until the present day. There 
is no dominant reservoir of metals at any epoch, although 
diffuse IGM gas holds a significant amount early on (mainly 
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because virtually all baryons are in this phase), while stars 
hold the most metals today. 

2.3 Missing Metals Problem 

Observational estimates have now been made for the metal 
content in various systems at z ~ 2.5, which as discussed 
earlier has led to a missing metals problem. Here we present 
a comparison of our model predictions with various obser- 
vations that try to account for the missing metals. 

At 2 ~ 2 - 3, FSB05 and BLP06 estimated that stars 
in Lyman break gala xies, particularly the BX sample of 
IShaplev et al. I (12005ft . contain about 18% of all metals. 
BLP06 go on to estimate a contribution of ~ 5% from so- 
called Distant Red Galaxies, though this may have some 
overlap with the BX sample jReddv et al. 1 ^005 ). Our pre- 
dicted metal fraction in stars is somewhat higher (« 30%), 
which may be explained if our typic al stellar mass o f LBGs 
is higher by ~ 50% than derived bv lShaplev et al. 1 {2005). 
This would easily fall within the substantial uncertain- 
ties of spectral energy distribution fitting techniques by 
which stellar masses are estimated. In general, simulations 
tend t o predict LBGs that have substantial older popula- 
tions llFinlator et al. 112006ft . and tend to be more massive 
than burst-dominated models would predict. Furthermore, 
there are completeness issues as our simulations include all 
stars in galaxies down to M* ~ 10 9 Mq, while typical LBG 
surveys are substantially incomplete at such masses. 

At z ~ 0, an observational census of metals by IPagefl 

(2002) finds ~ 50% in stars today depending on assump- 
tions about the IMF and stellar evolution, which again is 
in general agreement with our prediction of 43%. If we 
fix a stellar baryon fraction of 7% as our mod els predict, 
then lFinoguenov. Burkert fc Boehringer 1 (120031) determines 
Qz,* ~ 3 x 10 -5 , whereas we predict 2.5 x 10 -5 , again within 
uncertainties. In all, while more careful comparisons are war- 
ranted, the preliminary indication is that our stellar metal 
fraction agrees broadly with data, and shows that stars con- 
tain a substantial but not dominant fraction of metals in the 
universe at both high and low redshifts. 

Our simulations further predict that about 15% of met- 
als reside in cold star- forming gas in galaxies at z = 2.5. 
It has been suggested that Damped Lyman Alpha ab- 
sorbers (DLAs) may trace this cold dense gas. Yet de- 
spite the fact that such systems dominate the H I content 
of the Universe, t hey are seen to contai n just a few per- 
cent of all metals dProchaska et aL 2006). However, DLAs 
may .just be the tip of the iceberg nprochaska et al. I (12006ft 
and iKulkarni et al. I J2006ft find sub-DLAs (10 19 < N H i < 
10 20 ' 3 cm -2 systems) with above-solar metallicities, and ex- 
trapolate that such systems may represent 15 — 25% or more 
of the total metal budget. However, there is a bias towards 
selecting highly enriched sub-DLAs, so the extrapolations 
are highly uncertain. We will address sub-DLA metallici- 
ties in the next section. Another reservoir of star-forming 
gas is sub-millimeter galaxies, which are rare but massive 
and highl y enriched; from their gas c ontent and metallic- 
itv lBouche. Lehnert fe Peroux I (I2006af) estimated that their 
ISM may contain ~ 5% of the total metals. In general, it 
seems that the metals in star-forming gas are a relatively 
small fraction of the total metals, broadly agreeing with our 
predictions, although if sub-DLAs make as much of a con- 



tribution as claimed then perhaps the agreement is not so 
good. 

Combining these constraints, BLP06 determined that 
approximately 30% of metals are currently observed in 
galaxies (stars+ISM) at z = 2.2, which could plausibly in- 
crease to 50% with incompleteness corrections. This is quite 
consistent with our prediction of ~ 40% of metals being 
contained in galaxies with M* > 1O 9 M0 at these epochs. In 
particular, we do not predict, and our models do not require, 
a substantial amount of metals to be in some hidden form 
within galaxies. 

The remaining 60% of metals are not in the inner re- 
gions of galaxies, so have presumably been ejected through 
galactic outflows (note that the amount actually ejected is 
much larger, but many metals get reaccreted into halos). 
Of this, we predict that more than half of that is in the 
diffuse IGM at z ~ 3. This is substantially larger than pre- 
vious estimates. Using two-phase analytic modeling of the 
IGM, FSB05 estimate that less tha n 10% of meta ls are al- 
lowed in this diffuse phase (see also iPettini 1120041 and ref- 
erences therein). The discrepancy arises because the ioniza- 
tion corrections needed to obtain the metallicity are typi- 
cally done assuming the metals lie in gas having the density- 
temperature relation of the Lya forest. However, our simu- 
lations predict that a substantial amount of C iv-absorbing 
gas is heated to somewhat higher temperatures (see Figure 
9 of OD06) . Therefore we can hide more C IV in the diffuse 
phase (defined here as T < 10 4 ' 5 K) because our ionization 
fractions are lower than given by the Lya equation of state 
(see Figure 10 of OD06). Note that we still obtain a good 
fit to the observed 6-parameter and Ciii/Civ pixel optical 
depth ratio data, both measures of Civ absorber tempera- 
tures, despite this mildly shock-heated absorbing gas. 

The remaining 30% of metals in the Universe are di- 
vided between shocked IGM gas and halo gas at z ~ 2.5. 
FSB05 suggest that up to 90% of all metals may be hidden 
in hot galactic halo gas; our simulations do not support this 
idea. The amount in halo gas is difficult to quantify at high 
redshifts, but an estimate of metals in clus ters at z ~ 2.5 by 
iFinoguenov. Burkert fc Boehringer I (1200 3) predicts around 
10%, rising modestly to 20 — 25% at z ~ 0. Our predicted 
value for halo gas also evolves slowly as well, and the num- 
bers are similar, although a more careful comparison is war- 
ranted in order to critically assess how our definition of halo 
gas relates to gas in clusters or protoclusters. 

The fraction of metals in shocked IGM gas, clo sely 
related to the missing baryons (ICen fc Ostriker1ll999ft or 
WHIM gas, is re latively constant. S uch gas may be traced by 
Ovi absorbers iTripp et al. 1 120061 . though because Ovi is 
not the dominant ionization state of this gas (as we will show 
in H3.21 , the uncertainti es in the ionization corrections are 
large llTripp et aLj|200l J). Estimates of the metal content o f 
O VI systems bv lFinoguenov. Burkert fc Boehringer I d2003ft 
also show little evolution in their metal fraction, though they 
estimate around one-third of all metals are in this phase as 
opposed to our predicted values of 15 — 20%. 

So let us return to the original question: Where are the 
missing metals? According to BLP06 and our models, at 
least 50% are not in galaxies. The basic answer from our 
simulations is that they are mostly in gas that has been 
mildly heated by a combination of gravitational and feed- 
back effects. The enriched diffuse IGM is somewhat warmer 
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than expected from photoionization alone, the WHIM at 
these redshifts contains a non-trivial amount, and galactic 
halos also contain hot gas that has a non-negligible metal 
content. Hence metals are hidden in hot gas, but that gas 
is not confined to halos but rather spread among a wide 
range of overdensities (as we will show in ii3.ll . Overall, the 
simulations support a view where metals are rather evenly 
distributed among the various baryonic phases, which seems 
to be what observations suggest as well. 

2.4 Metallicity Evolution 

While obtaining a census of total metals is difficult because 
it requires detecting all metals in a given phase or extrapo- 
lating from available data, determining the mean metallic- 
ity is somewhat easier as it may have a characteristic value 
for a given phase of baryons. Hence the observations are 
somewhat better constrained, and in turn provide interest- 
ing constraints on our global enrichment model. 

The bottom panel of Figure shows the evolution of 
the mean mass- weighted metallicity in various phases, along 
with the global mean shown as the thick cyan line. The 
global mean metallicity tracks the diffuse IGM metallicity 
quite well down to z ~ 1, albeit higher by around 0.2 dex. 
This is largely because the diffuse phase contains most of 
the baryons in the Universe until z ~ 1, where WHIM gas 
begins to dominate. The global metallicity is around 2-3% 
solar at z = 2.5 and reaches a present-day value of almost 
one-tenth solar. 

It has been claimed that clusters represent a fair sample 
of the Universe in terms of galaxy formation processes, be- 
cause they have a stellar b aryon fraction similar to the global 
value fe.g. iRenzini lli"99g|) . However o ur models do not sup- 
port this idea (and neither do those of ICalura fc Matteucci I 
12004^1 . Our predicted hot halo gas metallicity is well above 
the mean global metallicity at all epochs, and does not 
evolve similarly. This is primarily because the interplay of 
outflows with cluster potentials allows larger halos to retain 
and/or reaccrete more of their metals. In other words, galax- 
ies are significantly affected by their environment not only 
through structure formation processes, but also through out- 
flow processes. 

Stars gain a mean metallicity of > 1/10 so- 
lar very early on, and in conjunction with a mass- 
metall icity relation whose slope is establi shed at early 
times iDave. Finlator fc Qppenheimer 1120061) . it is not sur- 
prising to find massive high-re dshift galaxies with solar 
metallicity IShaplev et al. 1120041) . The mean stellar metal- 
licity rises steadily from around one-tenth solar at z = 
6 to around one-half solar today. This latter value can 
be compared to observations, though there are difficul- 
ties in determining stellar populat ion metallicities ow - 
ing to degenaracies with age (e.g jTra ger et al. I 12000^1 . 
iFinoguenov. Burkert fc Boehringer I 12003) found a mean 
stellar metallicity of 0.6 solar in galaxy groups, which agrees 
nicely with our predicted value, assuming that stars in 
galaxy groups are representative stars everywhere in terms 
of their metallicity. 

Star-forming gas has a mean metallicity that evolves 
similarly to that of stars, only about a factor of two higher. 
Observational measures of gas-phase metallicity in galaxies 
generally rely on emission line tracers that arise in star form- 
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Figure 2. Moan gas metallicity weighted by star formation rate 
(solid), Hi column density (dotted blue), and free- free emission 
(dashed red), as a function of redshift. The first represents the 
mean metallicity of stars forming at a given redshift, the second 
corresponds roughly to the mean metallicity in DLAs, and the 
final one corresponds roughly to the X-ray weighted metallicity 
of intracluster gas. 



ing regions. Hence a more observationally relevant question 
is: What is the mean metallicity of stars forming today? This 
can be calculated as the SFR-weighted metallicity. 

Figure (solid line) shows the SFR-weighted metallic- 
ity as a function of redshift. This weighting is most closely 
related to the H II region emission line metallicity used to de- 
termine, e.g., the mas s-metallicity relation in galaxies (e.g. 
iTremonti et al. 112004) . At early epochs, the SFR -weighted 
metallicity is quite similar to the mean metallicity in stars, 
but by z = it is roughly twice the stellar one. At the 
present epoch, our model predicts that stars are forming on 
average at solar metallicity, although the average star has a 
metallicity half that. 

The metallicity of DLAs provides another constraint 
on global metallicity evolution. Their Hi column density 
weighted metallicity is seen to evolve slowly from one-tenth 
solar at z ~ 3 to one-fifth solar at z ~ .5. with some 
evide nce for lower metallici ties at z > 3 iPerou x et al. I 
2006 ) . In contrast, models bvlPei Fall fc Hauser 1 119991) and 
ISomerville. Primack fc Faber I 1200 lh . which match DLA 
metallicities well at z > 2, predict nearly solar metallicity 
at z ~ 0.5 for DLAs. 

In order to compare our simulations to this data, we 
determine an Hi column density- weighted mean metallic- 
ity, shown as the dotted line in Figure |2] This is ob- 
taine d by calculating the Hi ionization fraction assum- 
ing a lHaardt fc Madau I 1200 ll) background as described in 
OD06, and multiplying the Hi density by p _1//3 to obtain a 
(scaled) co lumn density for we ighting. In broad agreement 
with data (Pe roux et al. Il2006l) . the evolution is fairly mild 
from z ~ 3 — » 0.5, albeit somewhat higher than observed 
values. In contrast with models mentioned above, our pre- 
dicted DLA metallicity actually drops from z = 1 — > 0, and 
is distinctly sub-solar by the present epoch. 

It is worth noting that our computation of the H I frac- 
tion assumes optically thin gas; this is clearly an incorrect 
assumption in the highly dense regions that may be responsi- 
ble for DLA absorption. If we make the extreme assumption 
that every star-forming gas particle is completely neutral 
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(even though it is likely that at least some portion of the ISM 
is ionized), then the Hi column density- weighted metallicity 
goes up by 0.2-0.4 dex at z < 2. Obviously a more careful 
analysis is required; we leave this for future work. 

The discrepancy between previous predictions and data 
has spurred investigations of sub-DLAs as being responsible 
for holding h ighly enriched H I t hat m ay ma ke up the short- 
fall. I ndeed, iProchaska et al.1 i200rJ) and iKulkarni et al. I 
( 2006) find highly enriched sub-DLAs, and suggest that they 
may hold substantial amounts of metals, as discussed in H2.3I 
However, while their metallicities are quite high, this is offset 
by their lower column densities, so it is unclear whether they 
raise the global H I metallicity by a large amount. From our 
models' perspective, there is at most a minor discrepancy 
to explain: we predict mild evolution and significantly sub- 
solar metallicities for H I gas even at low redshifts, consistent 
with observations, though the somewhat higher metallicities 
overall predicted by our models could certainly accomodate 
some contribution from highly enriched sub-DLAs. 

A tight constraint for the metallicity of hot gas at z ~ 
is provided by intracluster media, which show a relatively 
uniform metallicity of around one-third solar. Effectively, 
this is an X-ray emission weighted metallicity. Figure [5] 
(dashed line) shows the metallicity of gas with T > 10 6 K 
weighted by its free-free emission, i.e. p 2 T°- 5 . The value at 
z — is around one-fifth solar, which is slightly below that 
seen in clusters. Note however that our simulation volume 
does not contain any actual clusters (our largest halo has 
a virial temperature of around 2 keV), and also free-free 
emission is a subdominant contribution to the emission at 
T < 1 keV where metal lines are prevalent. So the mod- 
est disagreement, while intriguing, should be regarded as 
preliminary. It is still worth noting that we predict little 
evolution in the emission-weighted metallicity of hot gas to 
z ~ 1 (in qualitative agreement with data, e.g. lRosati et al.1 
120041) . but at higher redshifts this metallicity drops rapidly. 

Shocked IGM metallicity remains relatively constant 
between one-hundredth and one-fiftieth solar, and at z < 4 it 
is very similar to the diffuse phase metallicity. The metallic- 
ity at z ~ in this p hase may be tested using O VI absorbers 
at low redshift (e.g. iTripp et al. I F2OO6I. Comparisons with 
simulations suggest that a typical metallicity of one-tenth 
solar in O vi-absorbing g as at low-z would give rise to a lin e 
frequency as observed jCen et al. ll200lllChen et al. 112003^ . 
These studies pasted on a metal distribution to hydro sim- 
ulations post facto. Since our simulations produce a lower 
mean metallicity than inferred by such studies, it bears test- 
ing whether our simulations can reproduce the observed O VI 
frequency directly from the simulated metal distribution. We 
are doing this comparison now. 

In summary, our models show broad agreement with 
the observed metallicity evolution as traced by stars, DLAs, 
and intracluster gas. Clearly more careful comparisons are 
warranted, but even these comparisons illustrate the power 
of metal observations in constraining outflow models, and in 
turn how such models can be used to coherently assemble 
information from a variety of objects and environments. 



0.06 



0.04 



0.02 



z = U 
z=l 
z = 2 
z = 3 



j : 1 

:" I 
.[ 'I 



2 4 
log p/<p> 




Figure 3. Histogram of metal mass in gas as a function of over- 
density (top panel), and temperature (bottom) at Z = 3,2,1,0. 
The area under the histogram is proportional to the total metal 
mass at that redshift. Shaded regions show the 2 = histogram 
for T < 10 4 - 5 K (top) and for p/p < 100 (bottom). The vertical 
axis values are arbitrary, but are proportional to the total metals 
mass Qz- 



3 METALS IN PHASE SPACE 
3.1 Evolution in Phase Space 

The evolution of metals cosmic phase space, i.e. overden- 
sity versus temperature, gives insights into galaxy formation 
processes. Metals are driven to lower densities through out- 
flows, and to higher densities through gravitational collapse. 
They are heated by gravitational infall and outflow-induced 
shocks, and cooled via radiative processes and Hubble ex- 
pansion. The path of metals in phase space over time is a 
key prediction of hierarchical galaxy formation models with 
outflows. Here we describe some general trends of the path 
of cosmic metals in phase space. We focus here on metals in 
gaseous form, leaving aside the metals locked into stars. 

Figure |3] shows histograms of metal mass as a func- 
tion of overdensity (top panel) and temperature (bottom), 
at z = 3, 2, 1 and 0. The overdensity here is taken from the 
particle densities, hence it is smoothed on the scale of the 
SPH kernel, which goes from hundreds of kpc in underdense 
regions to sub-kpc scales in the densest regions; it is the es- 
sentially the minimum scale over which density fluctuations 
can be modeled with this SPH-based code. At all redshifts, 
metals show a characteristic bimodality in overdensity, with 
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the diffuse peak holding the majority of metals. At z — 3, 
metals are prominent in the diffuse IGM with p/p < 1, but 
are also quite evident in a rather narrow peak at p/p ~ 200 
that likely arises because outflows stall as they intercept a 
dense IGM just outside halos. At z — 2, the diffuse IGM in- 
creasingly dominates the gas metal budget, but the peak has 
shifted towards higher overdensities. The trend continues at 
z — 1, and by by z = the peak of metals has shifted to 
fairly sizeable overdensities, around p/p ~ 200, and metals 
are now seen in very dense gas. 

The main trend from this is that metals shift to more 
overdense regions with time, because large halos are able to 
retain more of their metals. Note that in our simulations 
this is not because outflows escape more easily from small 
halos. Recall that in our momentum-driven wind model, the 
outflow velocity scales with the velocity dispersion. Hence 
all halos, to first order, lose a fixed amount of their ejected 
material. However, since the mass loading factor scales in- 
versely with the velocity dispersion, smaller halos eject a 
greater fraction of their star-forming gas, causing them to 
lose more of their total metal mass. Another factor is that, 
as star formation shuts off in the more massive halos, the 
outflows also subside, and the earlier-forming massive halos 
have more time to reaccrete their metals. 

The bottom panel shows the histogram with respect to 
temperature. Again, a sort of bimodality exists, although it 
is not obvious until z < 2. Very few metals have T < 10 4 K; 
although the underdense photoionized IGM extends to lower 
temperatures, the metals are restricted to the denser re- 
gions, and have a relatively small filling factor (OD06). 
Note that the cutoff at high temperatures is an artifact of 
our small simulation volume, whose largest halo at z = 
has a virial temperature of « 2 keV. While the amount of 
metals increases at all temperatures, the most notable in- 
crease occurs at T > 10 6 K. This owes primarily to the 
growth of large enriched potential wells from surrounding 
enriched WHI M gas, which i s then unable to cool their gas 
efficiently dKeres et al. 1120051) . However, this does not imply 
that cluster metallicities should evolve significantly, because 
as seen in Figure this evolution is driven mainly by an 
increase in the baryon fraction in the halo gas phase, while 
its metallicity remains fairly constant. 

The bimodality in both the density and temperature 
histograms suggests that the peaks may be straightforwardly 
related. But this is not so. To illustrate this, the shaded re- 
gion in the top panel of Figure |3] shows the low-temperature 
(T < 10 4 ' 5 K) portion of the z — density histogram, while 
the shaded region in the bottom panel shows the low-density 
(p/p < 100) portion of the z — temperature historgram. 
This shows that the high-overdensity peak is comprised pri- 
marily of cooler gas, while the major peak at moderate 
overdensities is mostly shocked gas. In the temperature his- 
togram, the low density gas dominates at 10 4 ' 5 < T < 10 6 K 
(corresponding to WHIM gas), while at T > 10 6 K and 
T < 10 4 ' 5 K one finds higher-density gas. 

In summary, gas-phase metals move into hotter, denser 
regions in the Universe with time, owing to the interplay 
between early outflows, gravitational reaccretion, and heat- 
ing by structure formation. At z = 2, the peak of metals 
is around the cosmic mean density, in diffuse photoionized 
gas. By z = 0, the metals have recollapsed into regions more 
than 100 times the mean density, with many of the metals 
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Figure 4. Phase space diagrams showing the metal mass (con- 
tours) and regions of phase space where the ions shown are the 
dominant state of carbon (left panels) and oxygen (right), at 2 = 2 
(top panels) and z = (bottom). Exceptions are Civ, Ov, and 
Ovi, which are shown when their ionization fraction is above 8% 
(or else they would not show up at all). Ionization states are 
colour coded as indicated to the left of the plot. 

shock heated to T > 10 5 K. Tracing this evolution directly 
with observations would provide a stringent test for cosmic 
enrichment models. 

3.2 IGM Metal Tracers 

How does one observationally catalog metals in various bary- 
onic phases? As discussed in H2.2I metals in ISM gas are 
probably most straightforward to catalog as they are visible 
in emissio n lines, though significant calibration iss ues re- 
main (e.g. lEllison. Kewlev fc MaHen-Ornelas"ll2005l) . Stars 
are also fairly straightforward to count, although the age- 
metallicity degeneracy requires spectroscopy to overcome. 
Metals in hot gas may be seen in X-rays, which can in prin- 
ciple pick up a significant fraction of metals at z = (Fig- 
ure [3J , but this is limited to low redshift systems having 
T > 10 6 K where the Galactic foregrounds are not over- 
whelming. The remainder of metals are probably most effec- 
tively traced through absorption line spectroscopy against 
background sources such as quasars, with the most com- 
monly utilized tracers being ultraviolet lines C IV and O VI. 
In this section we examine how such metal ions can trace 
the evolution of metals in phase space. 

Figure 01 shows a phase space diagram for gas from our 
simulations, with contours showing metal mass at z = 2 
(upper panels) and z — (lower). The evolution of metals in 
phase space as described in the previous section is evident 
here, as metals move to hotter, denser gas from z — 2 — > 
0. At low redshifts some metals do penetrate well into the 
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voids, affording at feast in principie the ability to trace low- 
density photoionized 1GM gas with metaf absorption fines, 
although the ionization ievefs are quite high. 

The coioured bands show regions in phase space where 
a given ion is the dominant ionization state for that metaf. 
Left panefs show carbon, right paneis show oxygen, and the 
colouring goes from singiy-ionized (in the iower right area) 
to highiy ionized (in the upper ieft area). For visibiiity, C IV, 
Ov, and Ovi are aiways shown when they have ionization 
fractions above 8%; this is done because those species are 
never the dominant ionization state, and hence wouid oth- 
erwise not show up on the piot. 

The key point to note on this figure is the overlap region 
between a given ion and the metal mass contours. This cor- 
responds to the range of densities and temperatures that are 
best traced by that ion. Of course, it is often possibie to de- 
tect muftipie ions for a given system, and the detectibiiity in 
practice is influenced by other factors such as the oscilfator 
strength and the efficiency of detectors in a given waveiength 
regime. Hence this piot is oniy intended to be iifustrative of 
broad trends. Generalfy, the bands show a characteristic pat- 
tern on the plot, with a diagonal rise from fow temperatures 
and densities owing to photoionization, and then a plateau 
at a particular temperature corresponding to the collisional 
ionization maximum for that gi ven ion. Note that t his piot 
was computed usin g CLOUDY jFertand et af. iFflfflSft with a 
Haardt & Madau (2001) ionizing background and assuming 
ionization equilibrium. The photoionized region is somewhat 
sensitive to our assumed ionizing background, but the colli- 
sional ionization plateau is not. 

From this plot it is straightforward to see why C IV 
(1548,155 lA) is a popuiar absorption feature for trac- 
ing fGM metafs, being an easiiy-identifiabfe doubfet at 
observationalfy-accessibie wavefengths. C IV overfaps the en- 
riched region of phase space from p/p ~ 10 — 1000 with 
T < 10 5 K at z = 2. As OD06 emphasize, Civ becomes 
a progressively better tracer of bulk 1GM metafs to higher 
redshifts. Conversefy, to fower redshifts it shifts to higher 
overdensities and becomes a poorer tracer of 1GM metafs. 
Cvi (33.7A) is more optimai for tracing photoionized IGM 
metals and C v (40. 3A) nicely tracks the bulk of metafs (and 
gas) in the WHIM regime, but both lie in the observation- 
ally challenging soft X-ray regime. At z = 2, Cm (977A) 
tracks the metals in cool gas around galaxies, but by z — 
it is predicted to be uncommon because there are few metafs 
left in such dense regions, as it has mostiy been consumed 
in star formation or shock-heated to higher temperatures. 
C II (1335A) traces gas in the inner regions of gafactic haios, 
which iikeiy invoives mufti-phase gas cfumps for which it 
becomes difficult to interpret fine ratios. Hence until signifi- 
cant advances are made in soft X-ray spectrographs, carbon 
is in practice oniy effective for tracing IGM metals at z > 2. 

Oxygen has the advantage of being a highly abun- 
dant element, and if starburst ejectae from Type II super- 
novae are primarily responsibie for enriching the IGM, may 
be afpha-enhanced refative to carbon. Furthermore, a wide 
range of oxygen fines are observationaily accessible and fairiy 
strong. The most common line used to see IGM metals is 
Ovi (1032,1037A), which traces a wide range of environ- 
ments from photoionized gas around the cosmic mean up 
to WHIM gas at overdensities around a thousand. Indeed, 
disentangling the photoionized Ovi absorbers versus the 



collisionally ionized ones is challenging dTripp et al. Il200ll . 
2006). At z > 2.5 this line star ts to become overwhelmed 
by th e thickening Lya forest iSimcoe. Sargent fc Rauch I 
l2004h . Ovu (21.6A) and Ovm (19A) t race the majority 
of WHIM metafs (e.g. IChen et al. 112003ft . as well as cover- 
ing the bulk of the photoionized metafs at fower redshifts. 
However, being that these are X-ray fines, the photoion- 
ized regions have too fow densities to provide detectabie 
absorption at current sensitivities, a nd the warm-hot sys- 
tems are only seen with great effort jNicastro et at. ||2005|. 
thoug h thi s detection is controver siat; see jKaastraet" akl 
12009) and iRasmussen et al. 1 12006)). Lower ionization lines 
like Ov (630A) and Oiv (787A) are interesting; although 
they do not dominate over a wide swath of phase space, it 
may be possible to detect them with UV instruments in con- 
junction with Ovi at moderate redshifts (z > 1), enabling 
detailed studies of the physical conditions and metallicity of 
T ~ 10 5 K absorbers. O ill (702A) traces denser, colder gas 
which at z > 2 contains substantial amounts of metafs, and 
11, like C 11, is limited to the regions immediately around 
star- forming gas. 

In the future it may be possibie to detect IGM metal 
lines in emission, though this will require a new UV- 
capable inst rument such as the p roposed Baryonic Struc- 
ture Probe llSembach et al. ( [20051 . The most natural ions 
to target from a practical standpoint are as usual O VI and 
C IV, because UV detectors are efficient down to ~ 1000A. 
However, it would be highly desireable to have detector ca- 
pabilities extending down to ~ 700A, so that a suite of 
O iii-O VI can be seen in relatively nearby large-scale struc- 
ture where the emission is not hopelessly dimmed by surface 
brightness effects. There is aiso devefopment underway for 
X-ray imaging with the D iffuse Intergalactic Oxygen Sur- 
veyor dOhashi et al. l f2006L Imaging offers the opportunity 
in principie to directiy trace the filamentary topology of 
WHIM gas, providing another interesting test of modefs. 

In summary, performing a full inventory of metafs in 
intergaiactic gas will require next-generation UV and X- 
ray spectrographs capable of detecting species that sample 
a wide range of ionization conditions. To this end, oxygen 
seems more promising than carbon, as O ill to O VI are all ac- 
cessible in the UV, and O VII and O Vlll are relatively strong 
X-ray lines, with oxygen being a highly abundant element. 
With statistical information on at least several of these fines 
for a farge number of systems, it will become possible to ob- 
servationaily characterize the metal distribution (along with 
the density and temperature distribution) of IGM gas, and 
directly test model predictions. For the time being, we are 
mostly limited to studying single ions in a small number of 
systems, and inferring the total metal mass from indirectiy 
determined ionization conditions or by comparing with mod- 
els such as the ones presented here. 



4 CONCLUSIONS 

Using cosmological hydrodynamic simulations of galaxy for- 
mation including galactic outflows, we have studied the dis- 
tribution and evolution of metals in various baryonic phases 
in the Universe. Our main conclusions are summarized as 
follows: 

• Metals are present in substantial amounts in a wide 
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range of environments from the earlier epochs studied (z ~ 
6) until the present. Diffuse 1GM gas dominates the metal 
budget early on, while stars and shocked gas hold most of 
the metals today. 

• The missing metals (outside of galaxies) at z ~ 2.5 are 
mainly in intergalactic gas that is mildly hotter than pure 
photoionization temperatures, having 10 4 < T < 10 5 K. 
However, this gas is not confined to halos, but is spread 
across a wide range of overdensities, with a peak around the 
cosmic mean density. 

• Comparisons of our simulation to observations broadly 
show agreement in the amount of metals in various baryonic 
phases, and show that metals are spread relatively evenly 
across all phases. 

• The global mean metallicity tracks about 50% higher 
than the diffuse IGM metallicity from z = 6 — > 1, after 
which it continues increasing to about one-tenth solar at 
the present epoch. Hot halo gas, which is predicted to have 
a luminosity-weighted metallicity of around one-fifth solar 
today, does not track the global metallicity, suggesting that 
clusters are unlikely to be representative samples of the uni- 
verse at large. This conclusions comes with a significant 
caveat, however, as our simulation volume is too small to 
produce actual clusters, so improved simulations are needed 
to confirm this. 

• Stars and star-forming gas quickly obtain a mean metal- 
licity of > O.IZq solar at z ~ 6, indicating early enrich- 
ment in galaxies. Today stars have a mean metallicity of 
one-half solar, while stars just forming now have typically 
solar metallicity. 

• Our models predict a global H I column density- 
weighted metallicity that increases slowly with redshift from 
z = 6 — > 1, and then actually decreases to z = 0, never 
exceeding around one-third solar. This is broadly consistent 
with observations of DLA metallicities, and does not require 
a substantial contribution from sub-DLAs. 

• Metals generally move from lower temperature and 
more diffuse regions at early epochs towards hotter, more 
overdense regions at low redshifts. This occurs because early 
blowout is very efficient at enriching diffuse gas, while at 
lower redshift outflows subside and large potential wells 
reaccrete gas. 

• Oxygen, particularly Oiii-Ovii, provides the most ef- 
fective tracer of metals across a wide range of cosmic phase 
space at z < 2, but a full metal accounting in the IGM will 
probably require next-generation instruments. 

Obtaining a complete census of metals in the Universe 
will require a concerted effort in studying a variety of sys- 
tems and environments at various epochs. Still, there is a 
realistic possibility of providing a full census of metals in 
the Universe in the foreseeable future. This is certainly a 
formidable challenge for upcoming facilities, but the payoff 
is large. By comparing such observations with cosmologi- 
cal simulations of outflows that track the interplay between 
various baryonic phases in mass, metals, and energy, it will 
be possible to stringently constrain models of galaxy for- 
mation and feedback. This will dramatically advance our 
understanding of how galaxies and intergalactic gas interact 
in the cycle of accretion and feedback, and provide for the 
first time a complete view of how baryons evolve in cosmic 
phase space from the dark ages until the present. 
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